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Gauge-Invariant QCD — an Approach to Long-Range 
Color-Confining Forces. 
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Department of Physics, University of Connecticut, Storrs, CT 06269 

We discuss the transformation of the QCD temporal-gauge Hamiltonian to a rep¬ 
resentation in which it can be expressed as a functional of gauge-invariant quark 
and gluon fields. We show how this objective can be realized by implementing the 
non-Abelian Gauss’s law, and by using the mathematical apparatus developed for 
that purpose to also construct gauge-invariant quark and gluon fields. We demon¬ 
strate that, in the transformed QCD Hamiltonian, the interactions of pure-gauge 
components of the gauge field with color-current densities are replaced by nonlocal 
interactions connecting quark color-charge densities to each other and to ‘glue’- 
color. We discuss the nonperturbative evaluation of these nonlocal interactions, 
which are non-Abelian analogs of the Coulomb interaction in QED, and we explore 
their implications for QCD in the low-energy regime. 


1 Introduction 

The problem that my student^] and I addressed in the investigation I am dis¬ 
cussing in this report, deals with the role of the non-Abelian Gauss’s law in 
generating a nonlocal interaction between gauge-invariant color-bearing parti¬ 
cle states — in particular, quark and multi-quark states. Our conjecture that 
this nonlocal interaction might play an important, perhaps a decisive, role in 
color confinement, is based on an analogy with Abelian gauge theories such 
as QED. In covariant-gauge QED, the spinor (electron) field if is a gauge- 
dependent quantity. It interacts with ‘pure gauge’ components of the gauge 
field A M , and that interaction is responsible for generating — perturbatively, 
order by order — the parts of the S-matrix elements that correspond to the 
Coulomb interaction between charges. On the other hand, when covariant- 
gauge QED is transformed to a representation in which Gauss’s law is im¬ 
plemented and if becomes a gauge-invariant charged field, the same nonlocal 
Coulomb interaction that is seen in the Coulomb gauge appears explicitly in 
that case, even though the covariant-gauge condition continues to apply.lH Quite 
generally, when only gauge-invariant fields are used in constructing the QED 
Hamiltonian, the interactions between charged fields and pure-gauge compo¬ 
nents of gauge fields vanish, and a nonlocal interaction between charge densi¬ 
ties — the Coulomb interaction — appears in their stead. In gauge-invariant 
QED, the only interaction besides the Coulomb interaction is between the 
transverse (gauge-invariant) part of the gauge field and the current density; 
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this interaction takes the form — fji(r)An (r)cZr, where An is the transverse 
(gauge-invariant) part of the gauge field, and j(r) = (r)dip(r). Since the 
current density has a v/c dependence in the nonrelativistic limit, the Coulomb 
interaction is, by far, the most important electrodynamic force in the low- 
energy regime. That fact provides a strong incentive for formulating QCD in 
terms of gauge-invariant fields, since then, in QCD too, the interactions be¬ 
tween the ‘pure gauge’ components of the gauge field and the color-current 
density must necessarily be eliminated and replaced by a nonlocal interaction 
between color-charge densities. We have obtained such a nonlocal interaction 
and, in the remainder of this paper, will discuss how it comes about and how 
it can facilitate our study of the low-energy regime of QCD. 


2 Implementing the Non-Abelian Gauss’s Law 


The strategy of using only gauge-invariant fields in formulating gauge theories, 
in QCD as well as in QED, faces major technical hurdles. The most daunting is 
the implementation of the non-Abelian Gauss’s law. The Gauss’s law operator 
for QCD is 

e“(r)=5 i n“(r)+ 5 r bc A t 6 (r),n i c (r)+j 0 “(r), (1) 


where 11“ is the negative chromoelectric field and, also, the momentum conju¬ 
gate to A“(r); j q is the color-charge density jg (r) = g ipH r ) V , and f abc 
represents the SU(3) structure constants. The Abelian Gauss’s law operator 
for QED is 

G(r) = <9*11, (r) + j Q (r) , (2) 


where j 0 is the electric charge density jo(r) = g ip^(r)ip(r). In QED, G(r) 
and <9,;IIj(r) are unitarily equivalent, which makes it far easier to implement 
Gauss’s law in QED than in QCD. We note, for example, that djllj(r) |£) = 0, 
or an equivalent equation, is easy to solve. Moreover, it is relatively straight¬ 
forward ±o explicitly construct a unitary transformation that relates G and 
Silh(r).l3 That unitary transformation enables us to transform solutions of 
dilb(r) |£) = 0 into solutions of G(r) |£) = 0, thereby implementing Gauss’s 
law. In QCD, however, Q a and 9,11“ cannot be unitarily equivalent, be¬ 
cause their commutator algebras differ. Whereas [<9jII“(r), 9;II*(r')] = 0, 


C) a (r),G b (r') = igf abc Q c (r)S(r - r') . □ Although 9jII“(r) |£) = 0 is just as 
easy to solve as its Abelian counterpart, we cannot make use of that solution 
to implement the non-Abelian Gauss’s law 


£“(r)|$) = 0. 


(3) 
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That fact requires us to use a fundamentally different strategy for implementing 
Gauss’s law in QCD. 

We have constructed states that satisfy Gauss’s law in QCDpand Yang- 
Mills theory, by solving an operator-valued differential equation, cm In that 
work, we defined an operator-valued quantity AJ{r), which we will call the 
resolvent gauge field , and which obeys the operator differential equation 


* J dr' [W(r), AJ (r 7 ) ] V? (r 7 ) + igf a ^ d A^ (r) J dr'[H d (r), A]( 
= -gf^ d At{ r)V d (r) 

+ E S ^r l3C f(nT A iV & (Mfjr)d,V/( r)) 


-EE 


(_l)*-y+* vKt ?$ t+1) /g A /“f n^ir) Mf v) ( r) diV?(T) 


-gf a ^( r) 


V" V'V_1 \ t rt t + r l b (’7) fP dx f 3x 7 di 
/ ^ ' 9 rj\(t+l)\ J (t) J{ri) d 7 

77=0 t=l 


«(r)^)(r)^W/(r)) 


(4) 


where 7?.“^ (r) is the product of functions of the longitudinal gauge field 


= fl *“ W ( r ) with A“(r) = [fM“(r)] . (5) 

m= 1 

Similarly, M^(r) is a functional of the resolvent gauge field 

v _ 

Mf v) (r) = 37“M( r ) with 3^(r) = r). (6) 

m= 1 

Eq.(|]) also contains the chain of SU(3) structure constants 

7 ) ( 7 ) 

the Bernoulli numbers B(g), and V~ ! (r). which represents an arbitrary vector 
field in the adjoint representation of SU(3). As shown in our earlier work, i 
Eq.(Q) is a reformulation of the non-Abelian Gauss’s law as a requirement on 
the resolvent gauge field Af{ r). In order to implement Gauss’s law, it is also 
necessary to -find an explicit solution of Eq. We provided such a solution 
in our work, El but, for lack of space, cannot repeat it here. 


yQ/37 _ j?a[l]f3b[l] ^b[l]a[2]6[2] yb[2]a[3]b[3] _ _ _ j‘b[ri—2]a[ri—l]b[ri—l] jb[ri 
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3 Construction of the Gauge-Invariant Fields 

The gauge field, and the resolvent gauge field, play a crucial role in determining 
the gauge-invariant spinor (quark) and gauge (gluon) fields. As we pointed out 
in earlier work,0 the Gauss’s law operator Q a { r), given in Eq. (0), and the ‘pure 
glue’ Gauss’s law operator Q a ( r) = dillf (r) + gf abc A b ( r)II?(r), are unitarily 
equivalent, so that Q a may be taken to represent <?“(r) in a different, unitarily 
equivalent representation. We refer to the representation in which Q a ( r) is the 
Gauss’s law operator, the C representation; and the representation in which 
0 represents the entire Gauss’s law operator, with the color-charge density 
Jq(i") included (though implicitly only), the Af representation. The unitary 
equivalence is expressed as El 

<3 a (v)=UcG a (v)Uc 1 , (8) 

where Uc is a functional of the gauge field, and the resolvent gauge field, as 
well as the color-charge density: lie = e c °e c with Cq and C given by 

C 0 = i J dr T“(r)jo (r) , and C = i j dr y a (r) j$ (r) . (9) 

Since the quark field ip trivially commutes with Q a ( r), ip is manifestly gauge- 
invariant in the AT representation. The unitary operator lie transforms the 
quark field ip so that the gauge-invariant spinor field in the C representation is 

V>ci(r) =Ucip{r)lic 1 = v c(r)i>{r) (10) 

where 

Fc(r) = exp(-i ff 3^(r)^) exp(-* 5 A a (r)^) . (11) 

and X h represents the Gell-Mann SU(3) matrices. And, since Vc(r) is a unitary 
operator, we can also readily see that the corresponding gauge-invariant gauge 
(gluon) field A b a (r) is given by 

[4uM £ ] = Vfc(r) [A b (r) £ ] V^r) + £ V c (r) d z V c ~\r) , (12) 

or, equivalently, 

4„(r)=4,(r) + fe-fg(r). (13) 

We can expand V’Gi( r ) and A^ u (r) to arbitrary order in g, and, in that way, 
we have verified that our gauge-invariant quark and gluon fields agree with the 
perturbative calculations of Lavelle, McMullan, et. a/fi-fo the highest order to 
which their perturbative calculations were available. 013 
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4 The Gauge-Invariant QCD Hamiltonian 

We can construct the gauge-invariant QCD Hamiltonian by systematically 
transforming the temporal-gauge QCD Ihimiltonian from the conventional C 
representation to the A F representation. E3 In that representation, ip(r) rep¬ 
resents the gauge-invariant quark field, and jg(r) the gauge-invariant color- 
charge density. In the Af representation, jft( r) therefore implicitly includes 
‘glue-color’ as well as the color of the bare quarks. 

The QCD Hamiltonian that results from the transformation to the Af 
representation is 

H= j dr [ in“(r)n“(r) + if 1 “-(r)I^“ (r)+^ t (r)(/3TO-*a i a i )V’(j')] + H' ■ 

(14) 

H' describes interactions involving the gauge-invariant quark field. The parts 
of H' relevant to the dynamics of quarks and gluons can be expressed as 

H' = Hj- A + H lr . (15) 

Hj-A describes the interaction of the gauge-invariant gauge field with the 
gauge-invariant quark color-current density, and is given by 

Hj-A = ~9 J drV’ t (r)a i ^V’(r)A£| i (r); ( 16 ) 


As in the case of QED, Hj-A couples the gauge-invariant gauge field to the 
current density j£( r) = gip^ and can be expected to have a v/c 

dependence for non-relativistic quarks, leaving the nonlocal Hlr as the dom¬ 
inant term in the low-energy regime. Hlr is the nonlocal interaction 


Hlr = H g -Q + Hq-q . 

A useful formulation of Hq-q can be given asi 

h Q-Q = \ J ^ r ^ x Jo(r).F ba (r,x)jo ( x ). 
where the Green’s function T ba (r, x) is represented as 


T ba (v, x) = 


47r|r — : 


- 2g f s ( 1 ) ba 


d y An 


47r|r-y| 




47r|y — x| 


(17) 

(18) 

+ 


- 1 - (_!)(" — !) 1 )g n . .fHn-i)S M a J 


dyi r^SV(y.) 


dy2 


47r|yi — y 2 | 


47T 


47r|r - yi|“ Gli 

——- 7 {y n ) : -y. (19) 

y(n-i)-y«| 47r|y„-x| 
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We make the following observations about J rba (v : x): The initial term resem¬ 
bles the Coulomb interaction. The infinite series of further terms consists of 
chains through which the interaction is transmitted from one color-charge den¬ 
sity to the other. Each chain contains a succession of ‘links’, which have the 
characteristic form 

= ( 20 ) 


The ‘links’ are coupled through summations over the S(„) indices and integra¬ 
tions over the spatial variables. All the quantities in Hq_q are gauge-invariant 
— the color-charge density j(j(r ) as well as the gauge field dg^(x). H g _Q — 


the other nonlocal interaction in Hlr — couples quark to gauge-invariant 
gluon color. In ti g -Q the quark color-charge density is coupled, through the 
same Green’s function J rba ( r, x), to a gauge-invariant expression describing 


‘glue’-color, given by K^(r) = gf dae Tr V c 1 (r)^Vc(r)^- 


^I,wn‘(r). 


A number of further observations about Eq. © have relevance for low- 
energy QCD dynamics. The first of these is based on a kind of ‘color multipole’ 
expansion of / J rho (r, x)j,° (x)dx — which appears as part of Hq_q as well as 
Hg-Q — about the point x = Xq: 


J dx jV ba (r,x 0 ) + XidiT ba (r,xo) + ^ Xi.Xj didjT^iy^o) + ■••jjo(x) 

( 21 ) 

where = (x — x $)* and di = d/dx.i. When we perform the integration in 
Eq. the first term contributes J rba ( r, x 0 ) Q a , where Q a = f dxjft(x) (the 

integrated ‘color charge’). Since the color charge is the generator of infinitesi¬ 
mal rotations in SU(3) space, it will annihilate any multiquark state vector in a 
singlet color configuration. Multiquark packets in a singlet color configuration 
therefore are immune to the initial term of the nonlocal Hq-q. Color-singlet 
configurations of quarks are only subject to the color multipole terms, which 
act as color analogs to the Van der Waals interaction. The scenario that this 
model suggests is that the leading term in Hq_q, namely Q b J 7ba (r Ql x Q )Q a 
for a quark color charge Q a at r 0 and another quark color charge Q b at x 0 , 
contributes to the confinement of quarks and packets of quarks that are not 
in color-singlet configurations. Moreover, assuming that T ba {r, x 0 ) varies only 
gradually within a volume occupied by quark packets, the effect of the higher 
order color multipole forces on a packet of quarks in a color-singlet configura¬ 
tion becomes more significant as the packet increases in size. As small quark 
packets move through gluonic matter, they will experience only insignificant 
effects from the multipole contributions to Hq-q , since, as can be seen from 
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Eq. (21), the factors Xi , XiXj, ■ ■ •, X^y ■ -Xy n y keep the higher order multi¬ 
pole terms from making significant contributions to fdx T ha (r, x)j'q (x) when 
they are integrated over small packets of quarks. As the size of the quark 
packets increases, the regions over which the multipoles are integrated also in¬ 
creases, and the effect of the multipole interactions on the color-singlet packets 
can become larger. This dependence on packet size of the final-state interac¬ 
tions experienced by color-singlet states — i. e. the increasing importance of 
final-state interactions as color-singlet packets grow in size — is in qualitative 
agreement with the characterizations of color tran sp arency and color coherence 
given by Miller and by Jain, Pire and Ralston. E30 


Another observation pertains to the series given in Eq. (|l9|) . To fully 
achieve the purpose of this investigation, it is necessary to evaluate X ha (r, x) 
nonperturbatively. The n th order term of J rba { r,x) — i.e. the chain with n 
links — has such a regular structure, that an explicit form can easily be writ¬ 
ten, without requiring knowledge of the lower order terms in the series. A 
nonperturbative evaluation of lF ba (r,x) is therefore far more accessible here 
than in the ‘standard’ gauge-dependent formulation in the C representation, 
in which it is much more difficult to identify and isolate the interactions that 
make the dominant contributions in the low-energy regime. 


In order to carry out a nonperturbative evaluation of Hq_q , either the 
J- ba (r, x) series has to be formally summed, or a more indirect method has to 
be found to achieve that objective. As of now, we have not yet addressed the 
evaluation of the gauge-invariant gauge field, which is an important component 
of X ba (r, x). We have investigated an SU(2) model of QCD under the assump¬ 
tion that gluon correlations can be neglected, so that the expectation value 
(0|Aj~ u (r)|0), in a state that implements the Gauss’s law C? a (r)|0) = 0, can be 
substituted for the operator-valued Ag u (r) in Eq.(^9[). Furthermore, we made 
the ansatz that the transverse {0|Ag| i (r)|0) can be modeled by the manifestly 
transverse SU(2) ‘hedgehog’ function (A|~ u (r)) = In this work, we 

showed how the SU(2) algebra enabled us to reduce the evaluation of J rba ( r, x) 
to the solution of a sixth-order differential equation,□ which circumvented the 
need for a perturbative evaluation of T ba (r, x). 


5 Discussion 

Abelian and non-Abelian gauge theories resemble each other in an important 
respect: The common thread that connects them is that the use of gauge- 
invariant fields eliminates interactions of pure-gauge components of gauge fields 
with charge and current densities, and replaces those interactions with nonlocal 
interactions among charge densities — the Coulomb interaction in QED, and 
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H L r , given in Eq. ©> in QCD. The structure of the Green’s function F ba (r, x) 
— a series of chains consisting of coupled links, in which the n th order chain is a 
product of n links — suggests features closely associated with QCD: flux tubes, 
‘string’-like structures tying colored objects to each other, etc. Further work 
will be required to determine whether these suggestive analogies are supported 
by detailed calculations. 
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